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Abstract 

^" 

Q^' The flavour and glue structure of the light scalar mesons in QCD are 

probed by studying the couplings of the / = mesons (t(600) and /o(980) to 
QJ I the operators qq, agC^ and to two photons. The Roy dispersive representation 

for the vrvr amplitude to(s) is used to determine the pole positions as well as the 
residues in the complex plane. On the real axis, tg is constrained to solve the 
^ _ Roy equation together with elastic unitarity up to the KK threshold leading 

■ to an improved description of the /o (980). The problem of using a two-particle 

. threshold as a matching point is discussed. A simple relation is established 

between the coupling of a scalar meson to an operator js and the value of 
the related pion form-factor computed at the resonance pole. Pion scalar 
form-factors as well as two-photon partial-wave amplitudes are expressed as 
coupled-channel Omnes dispersive representations. Subtraction constants are 
constrained by chiral symmetry and experimental data. Comparison of our 
results for the qq couplings with earlier determinations of the analogous cou- 
^ , plings of the lightest 1 = 1 and 7 = 1/2 scalar mesons are compatible with an 

?H \ assignment of the a, n, ao(980), /o(980) into a nonet. Concerning the gluonic 

operator a^G^ we find a significant coupling to both the a and the /o(980). 



1 Introduction: 

Exotic hadrons in QCD remain poorly understood theoretically. The recent discov- 
eries of the X, Y, Z states [Ij, for instance, in the charmonium spectroscopy was 
rather unexpected. Many of the expected states, on the other hand, which are associ- 
ated with gluonic excitations like hybrids or glueballs have not been unambiguously 
identified. The O"*""*" glueball is the lightest stable particle in the QCD spectrum 



1 



in the limit where all quark masses are sent to infinity. In this situation, its mass 
has been computed rather accurately in quenched lattice QCD simulations [2] to be 
slightly smaller than two GeV. In the presence of finite quark masses, the proper- 
ties of the glueball should remain relatively undisturbed provided mq > 1 GeV. In 
the physical situation, however, three quarks are substantially lighter than 1 GeV. 
Unquenched lattice simulations have been performed but the results are somewhat 
contradictory. The simulations of ref. [5] obtain near to maximal mixing between 
glueball and qq states and find that unquenching leads to a strong lowering of the 
masses. A similar effect of unquenching was observed for the 7 = 1 scalar mesons 
by several groups (e.g. [SI [6]). This picture, however, is not confirmed by the recent 
results from ref. [1] based on unquenched simulations with Nf = 2 + 1 and larger 
statistics who find gluball states very similar to the quenched ones. 

A possible scenario, suggested from using Laplace sum rules [3, E] is that there 
could be two mesons below 2 GeV with large glueball overlap. One of these could be 
rather light and possibly identified with the (t(600). Phenomenological implications 
of this scenario have been discussed in some detail in ref. [TO] . 

The classification of the lowest lying experimentally observed scalar mesons into 
a fiavour nonet is also not a completely solved problem [TT]. It has been proposed, 
for instance, that the ao(980) and the /o(980) mesons could have a specific status as 
weakly bound KK molecules [12] . This model simply explains their near degeneracy 
and their proximity to the KK threshold. It also seems able to explain the values 
of the 27 partial widths [T3|. Alternatively, it has been pointed out a long time 
ago that the mass pattern of the nonet below 1 GeV can be understood assuming a 
tetraquark fiavour structure [H] (see also [15]). 

The peculiarity of a nonet composed of the a, k, ao(980) and /o(980) is most 
clearly formulated in terms of 't Hooft's large Nc limit of QCD [16]. The masses, 
for instance, strongly deviate from the ideal mixing pattern predicted in this limitl§. 
This implies that in discussing the hght scalars, effects which are sub- leading in l/Nc, 
such as meson loops, ought to be taken into account. Modellings of meson loops 
effects can be found in the classic papers [HI [19] . More recently, a model from which 
an explicit l/Nc dependence can be deduced has been proposed [20]. Investigations 
in the ADS/CFT modelling of large A^^^ QCD have also been performed |21j . 

Experiments on radiative decays of the meson have been proposed [22| in order 
to clarify the fiavour structure of the light scalars. Such experiments have been 
performed and are planned to continue (see [23j for a review). The simplest way, 
however, to quantify the various aspects of the structure of the scalar resonances 

""^references to more recent work which incorporate, in particular, more realistic modelling of 
instanton effects can be traced e.g. from ref. |^ 

^In principle, dual ideal mixing is possible |15) . The scalars must then be either tetraquarks, 
i.e. exotics, or else the mass squared of the u-meson must be a decreasing function of the strange 
quark mass [17] which is unphysical. 
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would be via their couplings to a set of simple operators. The glue content, for 
instance, is best probed from the coupling to the gluonic operator agG"^. Similarly, 
the qq content is probed by the couplings of the scalar mesons to quark-antiquark 
operators. Such couplings have been considered for the 1 = 1 and 7 = 1/2 scalars 
by Maltman ^24j who suggested that their values can also be used for properly 
identifying the nonet. A lattice QCD result for the coupling of / = 1 scalars to 
ud is presented in Studies of couplings to tetraquark operators have also been 
recently undertaken [26| [27] . 

The cr(600) resonance is very unstable and does not give rise to a usual Breit- 
Wigner behaviour in cross-sections. Its existence has been demonstrated only re- 
cently [28j by making a combined use of experimental data and theoretical proper- 
ties of the TTTT scattering amplitude, which can be encoded into the set of Roy [21] 
integral equations. On the real axis, where the additional constraint of unitarity 
applies, the Roy equations were known as a powerful tool for analyzing experimen- 
tal pion-pion scattering data [30l [31]. New high precision experimental data on 
low energy pion-pion scattering [32l [33l [Ml [351 ES] have spurred renewed interest in 
these equations [3Zl [2Sl [2Sl [101 [II]- In ref. [M], the Roy equations are treated as 
a boundary value problem and exact solutions have been searched for numerically 
below a matching point s/sa = 0.8 GeV. 

When applied to resonances, the Roy equations are used for computing the 
partial-wave amplitude for complex values of the energy. The masses and widths 
of the resonances may be identified from the poles of the amplitude on the second 
Riemann sheet. The domain of validity of the Roy equations, as displayed in ref. [25] 
allows one to discuss both the a and the /o(980). The same poles which appear in 
the elastic scattering amplitude can be shown to also appear in two-point correlation 
functions of scalar operators and also in tttt matrix elements of these operators. The 
poles also appear in scattering amplitudes with a pion pair in the final state like 
77 —i- TTTT. The residues of the poles are also determined and can be interpreted in 
terms of couplings between scalar resonances and operators. In the present work we 
consider, from this point of view, the couplings of the scalar 7 = mesons a and 
/o(980) to the gluonic operator agG^ and to the quark operators uu + dd and ss. 
We will update the results that can be obtained for these couplings using the Roy 
equations combined with low-energy constraints from chiral symmetry. We will also 
consider the couplings to two photons, which were discussed in a similar framework 
in ref. [l2] . In that case, chiral constraints can be used as well as recent experimental 
data from the Belle collaboration ^3l [H]. 

The plan of the paper is as follows. We begin in sec. [2] by constructing solutions 
to the Roy equations in a domain which extends up to the KK threshold. This 
domain covers most of the /o effect on the real axis. We find that a very simple 
generalisation of the parametrisations used in ref. p8] is adequate for approximating 
the solutions. In sec. [3] we use these solutions inside the Roy integral representations 



3 



to perform extrapolations to the complex energy plane. We determine the resonance 
poles and their associated residues (sec. [S]). These results are applied in sec. IHto 
the determination of the scalar mesons couplings to two photons. For this purpose, 
we use the coupled-channel dispersive Omnes representation for 77 — >■ tttt, KK 
and the chirally constrained fits performed in [15]. The couplings of the scalar 
mesons to operators are finally considered in sec. O A complex plane definition is 
proposed from which a simple relation is obtained between the couplings and pion 
scalar form-factors computed at the resonance pole positions. Evaluations are made 
possible in this case by using chiral constraints for the form-factors in combination 
with coupled-channel Omnes representations [46j. 



2 Roy equation solution for t^{s) up to the KK 
threshold 



In order to improve the determination of the /o(980) properties, we begin in this 
section by constraining the J = S-wave amplitude to(s) to satisfy the Roy equation 
up to the KK thresholcH. The Roy equation reads 



Ret°(s) 



1 /""^ 
+ -/ ds' 



Imt°(.') 



+ KQis\s, 



+hiit\{s')Ki{s\ s) + \Yiitl{s')K2{s\ s) + dl{s) 



where Oq, 



Oq are the / 



0, 2 S'-wave scattering lengths. Detailed expressions for 
the kernels Ka{s',s) and the driving term d^{s) can be found in ref. [3^. Eq. ([1]) 
is supplemented with the non-linear, unitarity relation involving the inelasticity 
parameter ?7o(s) 



2ta^{s)tl{s)\ 



(2) 



with o",r(s) = yl — 4m^/s. 
one in the region s < 16m^ 



The inelasticity parameter 77° is rigorously equal to 
Based on experimental indications, we use here the 
approximation t^qIs) = 1 up to the s = 4m^. Furthermore in eq. ([T]), we use for the 
P-wave Im {t\) as well as the 1 = 2 partial-wave Im (tg) inputs taken from ref. [2H] 
i.e. satisfying the coupled Roy equations below 0.8 GeV and taken from experiment 
above. Imaginary parts of higher partial-waves, which enter into the driving term 
d^ are also taken from experiment. 



■^We neglect isospin breaking and take ttlk = {mK+ + ™i<'o)/2 
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2.1 Multiplicity of the solutions: 

Taking the matching point as Sm = sk = 4m|;, the Roy equation ([1]) admits a 
family of solutions |171 SHI |19] rather than a unique one. We will assume that the 
phase-shift at the KK threshold satisfies 

71<5k <Y ^ = SoisK), (3) 

which implies |171 HH HQ] a two-parameter family of solution^. In other terms, we 
must impose two conditions in order to select a unique solution. As one condition, 
we can fix the value of the phase-shift at one energy, for instance the value of 

6a = S',{sa), yiI = 0.8GeV. (4) 

In order to define a second condition, we consider the singularity of the derivative 
of the phase-shift at the matching point. For a generic Roy solution, the divergence 
depends on the value of the phase-shift at the matching point in the following 
way HZmg 



s)"-', a^?^^-2. (5) 



TT 



In our case, the matching point coincides with a two-particle threshold, we expect 
the derivative of the phase-shift to exhibit a square-root singularity 



A{sk-s)--^. (6) 



This divergence is weaker than the generic matching point divergence (E]) provided 
the threshold phase-shift is not too large, 

6^o{sk) < 225° . (7) 

We will assume here that this condition is fulfilled. In this case, we can use as a 
second condition that the phase-shift behaves as in eq. ([6]) close to the KK threshold. 
It is not difficult to work out the explicit expression for the coefficient A of the 
square-root singularity in eq. (jH]). For this purpose, let us consider the unitarity 
relation for Imtg in the region of the KK threshold 

ImtO(.) = a^{s)K{s)\' + 9{s - SK)MsMis)\' (8) 



One assumes that the following set of inputs are given: the two scattering lengths Oq, Oq, 
the phase-shift Sq{s) above the matching point, the inelasticity function rjQ{s) and, finally, the 
imaginary parts of the partial- waves ImiQ(s) and Imi^>-^(s). 
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where Qois) is the partial-wave tttt — ?► KK amphtude with / = 0, J = 0. The 
principal value integration in the Roy equation ([T]) generates singularities associ- 
ated with discontinuities of the derivative of ImtQ(s'). Finite discontinuities lead to 
logarithmic divergences upon integration. The square-root divergence is generated 
from the function 6{s' — SK)crK{s')- Performing the integration analytically in the 
neighbourhood of the threshold one easily finds that 

2 COS 26Ky/SK 

Once a solution is found for a given value of 6a, we can compute the over the 
experimental data in the range [sa, sk] and then search for the value of 6a which 
minimises this x^- In practice, the value of the phase-shift at the KK threshold, 
6k should be constrained by the data on both sides of the matching point. We can 
thus constrain both parameters 6a and 6k by fitting the experimental data using 
Roy equation solutions. 



2.2 Numerical approximations to the solution 

([1]) and by e(s) the difference 



Let us denote by 7^ [to] the right-hand side of eq. 
between the left and right-hand sides 



e(s)=7^[tO](.)-RetO(s) . 



(10) 



We construct numerical approximations to the phase-shift in the range 4m^ < s < 
4m|' using a simple modification of the Schenk parametrisation [50] compatible with 
eq. © 



tan 6q{s) 



N 



■Stt - So a^{s^) + f3 



with = 4m^ and ct^(s) = \/sk/s — 1. This representation involves polynomial 
parameters at plus 2 parameters sq and (3. The last factor generates a square- 
root divergence in the derivative of 6q as expected from eq. ([6]). In principle, the 
parameter /3 could be determined as a function of the known ([9]) coefficient A in 
front of the divergence. In practice, we have left it as a free parameter, adjusted such 
as to help approximate the solution for s close to sk but not necessarily reproducing 
the exact limiting behaviour for s = sk- We have checked that the correct order of 
magnitude for A is reproduced. 
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The N + 2 parameters in eq. flTT]) are determined from a variational principle, 
by minimising the integral over the error function squared 



4mi 



ds' 



(12) 



while fixing the two values of 5q{sa) and 5q{sk)- An exact solution corresponds to 
e(s) vanishing identically in the whole range [4m^,4m|f] and therefore to X/j = 0. 
We used routines from the MINPACK library [51] to determine the parameters in 
eq. f|TT]) which minimise Xr- We increased the number of parameters up to ten. 
With ten parameters one achieves an accuracy |e(s)| ^ 5 10~^ below the matching 
point. The behaviour of the error function is illustrated in fig. [H The figure shows 
that e(s) is an oscillating function which has a number of zeros approximately equal 
to the number of parameters in eq. flTT]) . The figure also illustrates how the error 
function evolves upon increasing the number of parameters which is suggestive of a 
convergence towards an exact solution. The accuracy is comparable to that quoted 
in ref. [38] below their matching point sa- Above the KK threshold, e(s) increases 
rapidly becoming ~ 10~^. In this region, this is similar to the results quoted in 
refs. 
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Figure 1: Error function (see eq. ( f7I}|) ) corresponding to an approximation of the 
solution (see eq. i fll]) ) with 8 parameters (dashed hne) and 10 parameters (sohd hne) 
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2.3 Inputs above the matching point: 

The behaviour of the inelasticity function r7o(s) close to the KK threshold is ex- 
pected to have a strong influence on the properties of the /o(980) resonance. By 
definition, 779 (s) is equal to the modulus of the tttt — t- tttt partial-wave S'-matrix 
element. Unitarity of the S'-matrix, 

|SnP = l-5^|Si„|2 (13) 

implies that //q = l^nl can be determined experimentally either a) by measuring the 
cross-sections of the various open inelastic channels or b) by measuring the cross- 
section for elastic scattering. The observation (by method (b)) that inelasticity 
sets in rather sharply at the KK threshold suggests that the KK channel should 
dominate the inelasticity below the rjr] threshold. The tttt — )■ KK amplitude with 
J = J = has been measured in high-statistics experiments [521 E31 [S4j- We will 
use here the results of ref. [52] because the results of [531 El] have been argued to 
necessitate some rescaling [551 EHj- The tttt — > rjrj amplitude has been measured 
in ref. [5^. Some experimental information on the tttt — )■ Att inelastic amplitude 
is also available. We will rely on the discussion of ref. ^6] who argue that the 
TTTT — )■ 47r amplitude is small in magnitude below 1.4 GeV and can be modelled by 
contributions from the /o(1370) and the /o(1500) resonances. 

Fig.[2]shows the experimental determinations of 7]^ based on the elastic amplitude 
from refs [58] and [59] . The result of the ii'-matrix fit performed in ref. [58] is plotted 
(dashed curve), which is characterised by a rather deep dip near 1 GeV. We also 
show the central value of the fit based on the inelastic channels (solid curve). The 
dip, in that case, is much less pronounced. The inelastic determination of r^g is 
actually not inconsistent with the elastic determinations of refs. [58] 159] within the 
errors. It has a /N = 1.6 with the data of ref. [58] and a /N = 0.4 with the 
data of ref. [59] (which is smaller than one because of the very large errors). For 
the phase-shift 6q above the KK threshold, we use the determination of Hyams et 
al. [5B] . It is in good agreement with other analysis of the CERN-Munich experiment 
(e.g. [56]) or the analysis of the CERN-Munich- Cracow experiment [59] below 1.5 
GeV. The energy region above 1.5 GeV is suppressed in the Roy equation because 
of the two subtractions. 

2.4 Inputs below the matching point 

In the energy range [5^,47711-], in which we fit the two parameters Sa and Sk, we 
combine the sets of data from Hyams et al. [55] and the data from Kaminski et 
al. [59]. The former data have much smaller error bars but it is likely that this is 
only because Kaminski et al. [59j have estimated their errors in a more realistic way. 
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Figure 2: InelsLsticity function rjois). The data shown are determinations from the 
elastic amphtude from refs. and J59|/ . The dashed hne is the K-matrix fit from 
ref. j58| . The solid line is a determination of rj^ based on experimental information 
on the inelastic channels KK, rjrj and Att. 



This is suggested by comparing the phase-shifts resulting from different analysis of 
the CERN-Munich experiment (e.g. [601 EH], see also the review [61j for detailed 
comparisons and further experimental references). We have therefore appended a 
weight factor of 1/4 to the of the data of Hyams et al. in the combined x^- Fo^' 
the S'-wave scattering lengths, we take the numbers quoted in the latest NA48/2 
publication 

a° = 0.2196 ± 0.0028stat ± 0.0020syst (14) 
al = -0.0444 ± O.OOOTgtat ± O.OOOSgygt ± O.OOOSQj^prj 

The results of fitting the combined data sets as described above in the region 
[s^,4m|;] varying the two parameters Sa and Sk are presented in tabled! We show 
separately the result corresponding to the two different determinations of the inelas- 
ticity function. We also show = /N (with = 10 data points) corresponding 
to the data of Hyam^ et al. [58j and to the data of Kaminski et al. [59]. The table 

^We remark that while the seems large, half of its value comes from the single energy bin 
with E = 0.99 GeV. 



9 



oo 


Sa 


5k 


^1581 
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(a) 


(80.9 ± 1.4)° 




2.7 


1.9 


(b) 


(82.9 ± 1.7)° 


(200l?o)° 


2.2 


1.3 



Table 1: Results for the two phases 6a and 6k from fitting the experimental phase- 
shifts in the range 0.8 GeV < \/s < 2mK with Roy solution functions corresponding 
to two different central values of the inelaticity function (see fig. |^. On the first 
line TjQ is determined from a sum over inelastic channels (shallow-dip shape), on the 
second line rj^ is determined from the elastic channel (deep-dip shape). 



shows that a better is obtained upon using the inelasticity function from the elas- 
tic data (deep-dip shape). This reproduces the observation first made in the recent 
analysis of ref. [H]. In that work, a variant of the three coupled Roy equations (de- 
rived from once-subtracted dispersion relations) have been considered in their whole 
domain of validity, i.e. up to y/s = 1.1 GeV and required to be satisfied withing the 
errors of the data. Their analysis favours a value for the threshold phase 6k some- 
what larger than the results of table [T] while their result for 6a_ is compatible with 
ours. Fig. 13] displays the curves for the phase-shift 6q corresponding to the fit results 
of table [H The figure also shows the phase-shifts from the Berkeley experiment [62] 
which were not included in the fit. Numerical values of the parameters describing 
the Roy solution phase-shifts (see eq. ffTTl) ) are given in the appendix. 

3 Poles and residues of the cr(600) and /o(980) 

Resonances correspond to poles of the tttt — ?► tttt scattering amplitude to(s) on un- 
physical Riemann sheets. These poles are also present in form-factors and correlation 
functions which involve currents which can couple to a pion pair in the S'-wave. We 
will consider only the second Riemann sheet here and recall a few standard formulas 
which enable one to perform the continuations If]. These formulas can be expressed 
in terms of the amplitude to(s). 

Let us start with the continuation of the amplitude itself. Its right-hand cut 
is associated with unitarity relations and has successive thresholds in s: 4m^, 16m^, 
36m^, KK, . . . The second sheet is defined with respect to the discontinuity relation 
which holds between the first two thresholds 4m^ < s < 16m^. Using the property 
of real-analyticity t^i^z*) = tQ*{z) (which results from T-invariance), it can be written 
as 

t[!(s + ie) - t°(s - ie) = 2(t^(s - ie)t%s - ie)tl{s + ie) (15) 
^More general formulas, for a four sheets situation can be found e.g. in ref. [53] 
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Figure 3: / = 5*- wave tttt phase-shifts: for y/s < 2mK the two curves represent 
solutions of the Roy equation corresponding to two different determinations of the 
inelasticity function rj^ (see fig. [^. 



where we have introduced 



(t"(z) = y^Aml/z - 1 



(16) 



(which satisfies a'^{s — ie) 
sheet extension of is 



i(7,r(s)). From relation (fT5|) . one finds that the second 



.0,11/ 



1 - 2a''(z 



17) 



which is easily seen to verify the continuity relation t^''\s 



s + ie). 



The poles of the T-matrix can now be determined by searching for the zeros of 
the denominator in eq. f|T7|) : 5*0(2:) = 1 — 2a'^{z)tQ{z) (which is the partial- wave 
5*- matrix). The derivative of Sq{z) 



Sl(.s) . ^(1 



2a'(z)tl(z)) 



=ZS 



is needed in order to determine the residues. Numerical results for the pole positions 
and the S-matrix derivatives are presented in table [2l The central values in the 
table correspond to the Roy solution associated with the deep-dip shaped tiq. This 
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(MeV) 




{zs) (GeV-2) 


a(600) 


(442+^) + i (274+^) 


- (0.75+ 


°i°)+.(2.20l°-) 


/o(980) 


(996^^,) + ^ {2AtV) 


-(1.1 





Table 2: Positions of the complex poles and values of the corresponding derivatives 
of the S-matrix Sq from the Roy integral representation of tg and the real-axis Roy 
solution discussed in sec. [2 



choice gives a result for the a position very close to that of ref. The errors 

were determined by varying the most significant parameters in the Roy equation i.e. 
the two scattering lengths a^, a^, the two phase-shifts 6q{sa), ^oi^K) (see table [T]) 
and the parameters of the /2 meson which dominates the driving term. We have 
also included the result of varying between the two different determinations of the 
inelasticity in the form of asymmetric errors. For instance, using the shallow-dip 
inelasticity, the value of the sigma pole position is located at : ^/s^ = 436 + z278 
MeV and that of the /o is located at y/sj^ = 983 + i36 MeV. The errors on the a 
pole parameters quoted in tableware smaller than those in [2E]: this can be traced 
to the fact than the range of variation for the phase Sq{sa) as determined from the 
fit using Roy solutions is smaller than the one estimated in ref. [2B]- 

4 Scalar meson couplings to two photons 

4.1 77 TTTT on the real axis: 

Informations on the couplings of the light scalar mesons to two photons can be 
extracted from the amplitudes 77 — ^ tt^tt^, tt+tt^. This may be performed in a 
model independent way by making use of the analyticity and unitarity properties 
of the partial-wave amplitudes /ijAA'l'^) which, as a consequence, satisfy Omnes- 
type [M] dispersive representations. A representation of this kind for /ig (s) was 
reconsidered recently [U] which makes use of a two-channel extension of the Omnes 
approach [HSl ES] • It should be valid in a range of energies up to one GeV where 
it is a reasonably good approximation to retain just two channels (vrvr, KK) in 
the unitarity relation. This representation involves also the 77 — )■ KK isoscalar 
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partial- wave amplitude fcg (s) and has the following form 



(19) 



lO.Bom/ \ 
, 0,Born/ N 



+ — 



s 

TT 



+ n{s) X 

ds' 



,(0) 

'k 



5'(o)s2 



'(0) „2 



-1/ 



Im 



ds' 



Ami 



is'Yis' - s 



-Imfl 



-1/ 



The right-hand side of this equation involves the 2x2 Omnes matrix f2, which 
encodes the effects of the final-state interaction. Its matrix elements are deter- 
mined from the T-matrix by solving (numerically) the set of homogeneous coupled 
integral equations which arise from combining dispersion relations and two-channel 
unitarity 

1 /""^ ds' 

%(s) = - / {T\s')nsW)),, (20) 

with E(s) = diag((j^(s), crA'(s)). One assumes asymptotic conditions on Tjj(s) (i.e. 
that (s) goes to zero and that the sum of the eigen-phase shifts goes to 27r |67] ) 
which ensure that eqs. fl20l) have a unique solution once initial conditions are specified 



fi,,(0) 



(21) 



These asymptotic conditions are rather close from the experimental values at y/s ^ 2 
GeV. Eq. (1191) also involves contributions from the left-hand cut of the partial-waves 
which are associated with singularities of the cross-channel amplitude 77r — )■ 77?. 
The leading singularity arises from the charged pion pole which is exactly calculable 
and labelled h^Q^°l^{s') in eq. f|T9|) (this term also dominates the amplitude in the 
soft photon limit). Singularities associated with multi-pion cuts are described more 
phenomenologically (but with reasonable accuracy) through the light resonance con- 
tributions, labelled /io'++(s') in the above formula. Finally, eq. f lTI?]) involves four 
polynomial parameters. These have been introduced by writing over-subtracted dis- 
persion relations, such as to cutoff integral contributions from higher energy regions. 
The polynomial parameters have been determined in ref. [35] from a chirally con- 
strained fit0 of the experimental data from ref. [13] (charged pions) and ref. [H] 



The fit was performed in an energy range -y/s < 1.3 GeV. For 1 — 2 amplitudes and for 
J — 2 amplitudes, single channel Omnes representations were used. Chiral constraints arise upon 
matching the dispersive and the chiral two-loop representations [5H1 ISH] from the fact that the p"* 
and certain chiral coupling-constants are known. 
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(neutral pions). In the present work, we use the vrvr phase-shifts obtained by solving 
the Roy equation below the KK threshold in association with the deep-dip inelas- 
ticity as discussed in sec. [2l As compared to ref. [15], this leads to small differences 
in the 77 — )■ tttt amplitudes localised in the region of the /o(980) peak. The values 
of the fitted parameters and the polarisabilities are not modified. 



4.2 77 — > TTTT in the complex plane 

Once the polynomial parameters are determined, the integral representations f|T9|) fl20l 
allow one to compute the partial- wave amplitude for complex values of s. In 

order to compute the second sheet extension one considers the discontinuity between 
the first two thresholds which reads, 

/io,++(s + ze) -hl_^^{s-ie) = 

2a''{s - ie)t°ois - ie)hl^^{s + ie), (22) 

such that the second sheet extrapolation is 

hoWA^) = i_2^.(,).o..^ • (23) 



The quantity of interest here is the decay width of the scalar mesons into two 
photons. Following Pennington [70] , it can be defined by first identifying the residues 
of the amplitudes t^'^^z) and /io+^+(2;) in terms of coupling constants 



32ntl''\z) = ^i^, /i^'++(^) = ^"""^"^^ . (24) 



,2 



pole Zg 



pole Zg 



The couplings gs'yy, Qs-k-k are expected to be complex numbers (see below). One 
can formally define the decay width by taking the usual relation between a coupling 
constant and the corresponding decay width 

which yields the following numerical results for the two-photons widths of the scalar 
mesons 

r.(600H2, = (2.08 ±0.20 (keV) 

r/o(980H27 = (0.29 ±0.21 ^^^y) _ 

The separation of the errors reflect the structure of the Omnes representation f[T9|) . 
The first error is associated with varying the subtraction parameters in eq. 0191) . 
i.e. it essentially reflects the experimental errors in the two-photon cross-sections. 
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The second error is associated with the uncertainties in the Omnes matrix elements 
coming from the tttt phase-shifts and inelasticities. Fig. H] compares our value for the 
sigma width with results quoted in the recent literature [TOl [TH UM UM US ESI ESI EZ] 
(see also [78]) which are all based on the complex pole definition. Evaluations using 
a Breit-Wigner definition can yield a somewhat different result (e.g. [IS]). In the case 
of the /o(980), which is a rather narrow resonance, the two definitions should give 
reasonably compatible results. The central value which we find in f l26p is practically 
identical to the one quoted in the PDG [80] . 



Present work 

Hofe richt_er (2 011) 

Mennessier (201 1) 

Mao (2009) 
■ 

Menn essier_(2008) 

Be rnabeu (2 008) 
Oiler (2008) 

Pen nington (2 008) 
O iler (2008 ) 

Pennington (2006) 

0.5 1 1.5 2 2.5 3 3.5 4 4.5 
Via ^ 27) (KeV) 

Figure 4: Recent determinations of the a — )■ 27 width from experimental measure- 
ments of 77 271 cross-sections. 



Let us finally quote the corresponding central complex values of the coupling 
constants gs-y-y, dsnn (in GeV), 

g^^^ = (-0.31 + zO.60) 10-2, g,^^ = 1.12 + z4.63 

= ( 0.38 - z0.02) 10-2, = 0.23 + i2.79 . ^ ^ 

It is striking that these couplings can be far from being real. It is difficult to find a 
general physical interpretation for the phases of the couplings but it is instructive to 
consider the case of a narrow resonance, i.e. when Im^;^ is small. In this situation, 
the Breit-Wigner approximation describes the amplitude in the region of the zero 
and the corresponding pole of the resonance 

S'oiz) ^ Ssik)'^ (28) 



15 



where k is the vrvr momentum, k = z/A — m^, and Sb is a slowly varying function. 
Neglecting contributions which are quadratic in Imz^, this representation gives the 
derivative at z = zs as 

. ^ zexp(2^^g(Re..)) 
2Im zs 

Using this in the expressions for the residues, the coupling gsmr gets expressed as 

9snn _ -exp{-2i5^iRezs))lmzs 

327r a^{Rezs) ' ^ ' 

i.e. the phase of gsmr is given in terms of the phase-shift at the resonance mass 

gs.. = |^?5..|e<^-^o(^<^-^) (31) 

and vanishes only in the absence of any non-resonant background phase. In the case 
of the coupling gs-y-y one finds, at leading order in Imzs, 

gsirngs-r-y - 2iexp{-2i6Q{Rezs))hQ{Rezs)lmzs (32) 

i.e. (using (13T]) ) 

Phase [gsyy) = Phase (/io(Re zs)) - 5^{Re zg) (33) 

which vanishes modulo tt when Re zs is in the region of applicability of Watson's 
theorem. These narrow width estimates for the phases provide a reasonably good 
approximation for the /o(980) when its mass is located below the KK threshold 
(which is not the case for our central value). 

5 (j(600) and /o(980) couplings to gluon and quark 
operators 

5.1 Definitions 

The cr(600) and /o(980) mesons have the same quantum numbers J^'" = 0^~^ as the 
vacuum which are also those expected for the lightest glueball. One can characterise 
the gluon content of a scalar meson from its coupling to the gluonic operator a^G^. 
One may also consider the trace of the energy- momentum tensor operator, 6'^, which 
is proportional to OsG^ in the chiral limit. Correspondingly, two coupling constants 
C^^, Cg (with mass dimension) can be introduced 
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where S is either the a or the /o(980) meson. We will also consider matrix ele- 
ments associated with scalar quark-antiquark operators. It is convenient to use a 
normalisation which remains well defined in the chiral limit 



\uu 



+ dd\S) = V2B0 Q", {0\ss\S) = Bo Clf . (35) 



with Bq = — limm^^o {0\qq\0) / . With this convention, the couplings are renor- 
malisation group invariant in the chiral limit. 

At first, it is necessary to clarify the meaning of such matrix elements since scalar 
mesons are resonances and not stable one-particle states. One may use a complex 
plane definition, which is rather natural here as it applies equally well to broad 
resonances like the a or to ordinary narrow resonances. A simple relation between 
the couplings and pion scalar form-factors can be derived. For this purpose, let 
us consider two-point correlation functions 

%(s) = z J d^xe'P^{0\Tjs{x)jsmO) (36) 

where js{x) is one of the scalar operators considered above. The correlator Yljj 
satisfies a Kallen-Lehmann representation (e.g. [8T] ) 

(written here with three subtractions) in which the spectral function is given as a 
sum over a complete set of states 

(27r)^ J]5^(p„-g)|(0b5(0)|n)r = ^^(go)p,,(g^) • (38) 

n 

The discontinuity of Ujj across the real axis in the range 4m^ < s < 16m^ is 
generated by the two-pion states n = vr^vr" in the sum (1551) and it can be written as 

%(s + 2e) - %(s - ze) = ^a-{s - ie)F,{s - ie)F,{s + le) . (39) 



Here, Fj is the form-factor associated with the two-pion matrix element of js 

(0|j5(0)|7rXp)vr^(p')) = 5^' F,{{p + p'f) . (40) 



In deriving eq. fl5I?l) one makes use of the fact that Fj{s) is itself a real-analytic 
function. It has a cut along the positive real axis, and its discontinuity in the range 
[4m^, 16m^] reads 

Fj{s + ie) - Fj{s - ie) = 2a^(s - ie)t°o{s - ie)Fj{s + ie). (41) 
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From the discontinuity relations ( 15^ ( HTj) . it is simple to deduce the second sheet 
extensions of the form-factor 



and that of the correlator 11^^ 



-{z) {F,{z)f 



njjWU-^^^^^^^^, (44) 



n-(^)-%(.) + - ^_— ^ . (43) 

These expressions show that the form factor and the correlation function on the 
second Riemann sheet have exactly the same poles zs as the T-matrix. Considering 
the residue of the pole provides a natural identification for the resonance couplings 
ms\S). 

m3s\S)f 
Zs- z 

which thus get expressed in terms of the vrvr form-factor evaluated at the position of 
the pole, 

One can verify that the interpretation of residues in terms of coupling constants 
satisfy consistency conditions. For instance, one expects the residue of the form- 
factor Fj\z) to involve the product of the two couplings (0|js'|S') and gsnn in the 
following way 



Ff(z)| , =m0l211^_ (46) 



VS{zs - z) 

It is easy to verify that this expression can be exactly recovered using formu- 
las f|T7|) . fH2|) ■ fH3|) for the second-sheet extensions together with the definition of 
gg^^ from the residue of tg'^^(z) and the definition of (OljsIS') from the residue of 

nll(.). 

In the limit of narrow resonances, one can express the couplings Cg in terms of 
the form-factor Fj evaluated on the real axis. For this purpose, one can write Fj in 
the neighbourhood of the resonance position as a function of the momentum k 

FA^) - (47) 

s 

displaying explicitly the pole on the second sheet. If the pole is close to the real axis 
we can expand the function (f)j{k), 

(Pjiks) = (pjiReks) + i{lmks)(j)'{Reks) + ■■■ (48) 
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which, to lowest order in Im ks leads to the approximation 

Fj{zs)^lFj{Rezs) . (49) 

Using also the expression for the derivative of the ^-matrix in the narrow width 
limit ( 129|) one can express the couplings in terms of quantities evaluated on the real 
axis 

mjs\S)f ^ ^cr^iMDMsTs {e-<^^s) f^^mI))" , (50) 

using Rezz — M|, Imzs = MsTs- This expression shows that the squares of the 
couplings Cg must be real numbers in the narrow width limit, provided Ms is in 
the region of applicability of Watson's theorem. The couplings themselves can be 
either real or pure imaginary depending on whether the phase shift and the phase 
of the form-factor are equal or differ by tt. 



5.2 Numerical results 

Analyticity and unitarity allows one to derive Omnes representations for the form- 
factors, analogous to those for the 77 — )■ tttt amplitude but much simpler because 
of the absence of a left-hand cut. Let us briefly recall the derivation. Let F{s) be a 
two-component vector formed from the pion and kaon form-factors, 

'm = {F^{s),^Fj'{s)) (51) 

and multiply it with the inverse of the Omnes matrij|f| 

G{s) = n-\s)F{s) . (52) 

This multiplication removes part of the right-hand cut i.e. the components of G{s) 
have a right-hand discontinuity which vanishes in the range 

ImG(s) ~ 0, Ami < s < S2 (53) 

where S2 is the point above which two-channel unitarity is no longer a good approx- 
imation. By construction, the components of 0(s) behave as 1/s when s — )■ 00 and 
a similar behaviour is expected from the from-factors, such that G{s) should satisfy 

^The determinant of the Omnes matrix can be expressed in analytical form: detr2(s) = 
{ilZl with <P{s') = e{Aml - s')S^oW) + e{s' - 4m|)5,,^^^(5') which shows 

that it does not vanish. 
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a once-subtracted dispersion relation. In terms of F, it reads 

F{s) = n{s) 



s 

+- 



ds' 



Im {n-\s')F{s')) 



(54) 



In the range s << S2, the energy dependence of the integral may be neglected and 
one ends up with the following representation for the form-factors 




a + a's 
+ 



(55) 



As the discussion above shows, it is valid for s << 82- Such representations were 
used and discussed in detail in ref. [IB]. In order to determine the polynomial 
parameters, one can rely on chiral symmetry PB]. As a first approximation, one 
can use the chiral expansions of the form factors at order and determine the 
polynomial coefficients by matching the O(p^) values of Ff{0), Fj^ (0) with those 



3s 


FJ{Q) 


F7(0) 


Ff{Q) 


^f(O) 


muUU + rriddd 


ml 





1 2 
^ml 





rrisSS 








2 12 





0^ 


2ml 


1 


2ml. 


1 



Table 3: Pion and kaon form-factors associated with various operators js- The table 
shows their values at s = and the values of their derivatives at leading chiral order. 



of the Omnes representation. These 0{p'^) values are recalled in table [31 The 
representation fl55|) then allows one to compute the form-factors for complex values 
of s (with |s| < S2) and thus determine the values of the couplings between scalar 
operators and scalar mesons from residue relations like psl) . 

The numerical values of the absolute values of couplings (the phases will be shown 
later) of the a and /o(980) mesons to the qq operators obtained in this manner are 
collected in table [H In this table, the first error refiects the infiuence of higher 
order chiral corrections in the polynomial parameters. We have estimated that the 
order of magnitude, relative to the OijP') values, should be ~ 30% for the corrections 
proportional to m^, and neglected the corrections proportional to m^^. As expected, 
a larger uncertainty is generated for the /o(980) than for the a. The second error 
is associated with the uncertainties in the tttt and KK T-matrix as refiected in the 
Omnes matrix elements. 
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ff(600) 


/o(980) 


1 /^uu 1 
\'^S \ 


(MeV) 


206 ± Atl 


82 ± 31^^^ 


\r^ss 1 
\'^S \ 


(MeV) 




146 ± 44+^^ 



Table 4: Absolute values (in MeV) of the couplings of the a and /o(980) mesons to 
scalar qq operators as defined in eq. (135)) . 



A previous estimate of the qq coupling of the a meson, using Breit-Wigner ap- 
proximations, was given in ref. [82] in the form {0\dd\a) = >/2/3Bo/x, with x = 20 
GeV~^, which is significantly smaller than our result. Some results for the couplings 
of the 1=1 and J = 1/2 mesons to qq operators can be found in the litterature. 
These resonances are reasonably narrow, such that various definitions should be 
equivalent and we can compare their values to those we found for the / = mesons. 
We normalise the couplings of these mesons in accordance with eq. fl35|) 

{0\us\K*) = B,Cl\, {Q\ud\ao) = BoC^^^ . (56) 

An evaluation of the ao(980) coupling was performed in ref. [21] using finite-energy 
sum rules (see also ref. [S3] ) . Converted to the normahsation of eq. , the result 
of [21] reads, 

|Cat(980)l = 197±37MeV (57) 

which is remarkably similar to the coupling of the a meson C"" in table HI The 
coupling C*"^'^(.ggg-) is related to the coupling introduced in ref. [H] by C^g^^g^Q) = 
Acm and can be estimated from its relation to the low-energy chiral coupling con- 
stants [84j, eventually supplemented with large or chiral sum rule constraints [S51 
[86] . These approaches yield values in the range (^"^''(ggo) = [120,200] MeV. An 
unquenched lattice QCD calculation has also been performed [25] which gives: 
^ao1;980) ~ [^^^y 340] MeV. These values should not be compared too litteraly to the 
preceeding ones because they correspond to unphysical pion masses / mPj^^'^ > 5 
and only two dynamical flavours. An estimate for the k meson coupling C"** can 
be made following a similar approach to that used here for the a meson. One can 
compute the position of the complex pole and the corresponding value of the 5- 
matrix derivative from the Roy-Steiner equations [87]. The central values which one 
obtains in this way are 

~ (658 + i 277) MeV, sl (zs) ~ (0.59 + i 2.03) GeV~^ (58) 

The coupling can then be defined in terms of the Kir scalar form-factor evaluated 
at Zs (see ref. [88], appendix C) and this gives 

|Q^800)l^l56MeV. (59) 
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Comparing now the couplings of the 1 = mesons from table H] to those of the 
1 = 1, 1/2 mesons one observes that the values of C^", C^^, C%{9so)^ ^ao{9so) 
rather similar, the relative differences do not exceed ~ 20%. This is compatible 
with an assignment of the mesons a, k, /o(980), ao(980) into a nonet. Results on 
the couplings of the heavier scalar mesons ao(1450) and -ft'Q(1430) are also available. 
Ref. [21] gives 

|Coa450) I = 284 ± 54 MeV, |Q^*(i43o) I = 370 ± 20 MeV . (60) 

The result for the ao(1450) was obtained from a finite-energy sum rule and the one 
for the i^'o(1430) from a one-channel Omnes representation. An evaluation using a 
two-channel representation and complex pole definition was made in ref. |SB] which 
gives |C*^**(i43o) I — 282 MeV. With the normalisations used here, the couplings of 
the ao(1430) and -^"0(1430) to quark-antiquark operators seem to be significantly 
larger than those of the light scalars. 





a(600) 


/o(980) 


\Cl\ (MeV) 
ICf^l (MeV) 


197± 15tf 
472 ± 


114±44tf 
227 ±41+^^ 



Table 5: Absolute values of the couplings of the a and /o(980) to the gluonic oper- 
ators 9^, and asG^ . 



Finally, one can compute the couplings of the light / = scalars to the energy- 
momentum trace operator 6'^ using the chiral results for the associated form-factor 
at s = from table [3] . The results are shown on the first line of table O One finds 
that both the a and the /o(980) display a significant coupling to the 6^ operator. 
The trace of the energy-momentum tensor has the following exact expression in 
QCD [93j with three heavy flavours integrated out 

0', = ^G^^c^^-^ + (1 + ^Ug)) E • (61) 

q=u,d,s 

This expression allows one to disentangle the a^G^ part from the qq one if one uses 
a perturbative approximation for the /3 function and for the anomalous dimension. 
The results shown in table O for Cg'^ correspond to a leading order approximation. 
Our results for C| may be compared with the Laplace sum rule evaluation [S] 

= [272, 329] MeV . (62) 
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However, one should keep in mind that in the calculation of [S], the spectral func- 
tion lmlljj{s) corresponding to the operator js = 0'^ is approximated by a simple 
delta function. Fig. [5] shows our result for this spectral function based on using 
two-channel unitarity and physical vrvr scattering inputs. It displays a peak corre- 
sponding to the /o(980) resonance, while the a resonance does not show up as a 
clear enhancement, but generates a broadening of the /o(980) peak at low energies. 
It is then plausible that the value ( 1^ should be compared with the sum + Cj^ 
from table O the agreement is then rather reasonable. 
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Figure 5: Spectral function of the Hjj correlator with js = 0'/^- The long-dashed and 
short-dashed curves are the contributions from the vrvr and KK intermediate states 
respectively. 



Finally, the central values of the phases of the couphngs Cg are shown in table O 
In the Breit-Wigner approximation, one expects the phases to be either zero of ±90° 
(see fl50l) ). The actual values are often not too different from this approximation. 



6 Conclusions 

We have considered several properties of the light scalar isoscalar mesons a and 
/o(980) using definitions which rely on the positions of the poles in the complex 
plane and their associated residues. This approach allows one to deal with a broad 
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a 


/o(980) 


uu + dd 


28.2° 


89.1° 


ss 


-80.2° 


-14.2° 


K 


87.2° 


-34.6° 





Table 6: Central values of the phases of the couplings Cg. 



resonance like the cr in a well defined way. In order to compute the positions of 
the poles and the residues,the Roy integral representation for the vrvr scattering 
amplitude was used. On the real axis, we have started from the Roy equation 
solutions of ref. |3H], which use a matching point a/s^ = 0.8 GeV and construct an 
extended solution which, for the 5- wave tg, has a higher matching point ^/s^ = 2mK 
such as to improve the theoretical constraints on the /o(980) meson properties. 
In order to constrain the value of the S-wave scattering phase-shift at the KK 
threshold and discriminate between different shapes of the inelasticity, corresponding 
to different experiments, we perform fits of the phase-shifts below the KK threshold 
based on the Roy solutions. We find that the solution corresponding to a deep- 
dip shaped inelasticity has a better than that corresponding to a shallow-dip 
shape. This is in agreement with the observations of ref. ^41j. The properties of the 
/o(980) resonance, as expected, are particularly sensitive to the central value of the 
inelasticity. The results based on this Roy representation of the amplitude for the 
second-sheet pole positions are in table El 

As a first application, we have re-determined the scalar to two photons cou- 
plings Qs-yy, following the methodology first advocated in ref. [12], and based on the 
determinations of the 77 — t- tttt amplitudes from the recent experimental measure- 
ments [121 SI]. The result found for the a is somewhat smaller than that originally 
given in ref. |12]. As a second application, the couplings of the a and /o(980) mesons 
to scalar operators, which can be formally denoted as {0\js{0)\a), (0|j5(0)|/o(980)) 
were defined and evaluated. Choosing js = (uu + dd) j -\/2, js = ss these matrix ele- 
ments provide a quantitative measure of the quark-antiquark contents of the scalar 
mesons, while choosing js = is a measure of the glue content. A simple, gen- 
eral relation can be established between such couplings and the value of the pion 
form-factor associated with the operator js computed at the position of the reso- 
nance pole, FJ'^{zs)- This relation is given in eq. f HSjl in the general case and in 
eq. fl50|) in the limiting case of a narrow resonance. Such form-factors are known 
to be calculable from a coupled- channel Omnes representation [4Sj which should be 
valid in a complex energy range which accomodates the a as well as the /o(980) 
resonances. The polynomial parameters in such representations are constrained by 
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chiral symmetry, for both the qq and 6'|^ operators, and can be estimated from the 
leading order chiral Lagrangian [IB]. In principle, matrix elements of other types of 
operators, for instance tetraquark operators, could be addressed in the same way. 
The values of -Fj(O) and Fj(0), in such not predicted from chiral symmetry 

but could be obtained e.g. from lattice QCD. 

The numerical results for the qq coupling constants of a and the /o(980) mesons 
are shown in table HI The couplings are not particularly suppressed but it would be 
interesting to compare them with couplings to tetraquark operators. The couplings 
can also be compared to the analogous couplings of the 1 = 1 and 7 = 1/2 mesons 
to the ud and us operators respectively for which estimates can be found in the 
litterature including one calculation in lattice QCD [25]. This comparison supports a 
nonet assignment of the a, k, ao(980), /o(980) mesons. Our results for the couplings 
to the gluonic operators and asG"^ indicate that both the a and /o(980) couple 
significantly to such operators as well. 
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Appendix 

We show below central values of the parameters of Roy solutions for the phase-shift 
5q(s), for a ten-parameter approximation according to eq. (ITTi) . corresponding to 
two different central values of the inelasticity function ?7q, see sec. [2l 





1]^: deep-dip 


t]q: shallow- dip 


So 


0.724237452 


0.736126142 


/3 


0.104114178 


0.360170063 


ai 


0.140785825 
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a2 
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-2.8956852410-6 
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